Evidence for the higher twists effects in diffractive DIS at HERA 
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We study a twist decomposition of diffractive structure functions in the diffractive 
deep inelastic scattering at HERA. At low Q 2 and at large energy the data exhibit 
a strong deviation from the twist- 2 NLO DGLAP description. It is found that this 
deviation in consistent with higher twist effects. We conclude that the DDIS at 
HERA provides the first, strong evidence of higher twist effects in DIS. 

I. INTRODUCTION 

The QCD description of the diffractive deep inelastic scattering processes ep —> epX 
(DDIS) is based on the series expansion of the scattering amplitudes in the inverse powers 
of a large scale Q 2 , denned as a negative squared four-momentum transfer from the electron 
to the proton carried by the virtual photon 7*. In the leading twist- 2 approximation the 
diffractive proton structure functions can De calculated using diffractive parton distri- 

bution functions (DPDFs) due to the Collins factorization theorem [l| , whereas the DPDFs 
dependence on the hard scale is governed by the celebrated DGLAP evolution equation. De- 
spite of great efficiency of this approximation in the data description this approach has an 
obvious limitation that follows from negligence of the higher twists contributions. Certainly, 
the higher twists contribute at any energy scale and become relevant for data description 
below some virtuality Q 2 , which depends on the process and required precision. In this 
presentation we point out that in the case of DDIS the DGLAP description breaks down at 
the scale Q 2 ~ 5 GeV 2 and to show that these deviations are consistent with a higher twists 
contribution. 
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Fig. 1. Left panel - kinematics of the DDIS scattering. Right panel - the x 2 / d.o.f. for NLO 
DGLAP and NLO DGLAP + HT fits to ZEUS LRG data [2] with Q 2 < Q 2 min . 



II. CROSS SECTION AND THE DGLAP DESCRIPTION 

The DDIS is an quasi-elastic electron-proton scattering process e(k)p(P) — > 
e(k')p(P')X(Px) in which the final hadronic state X with four-momentum Px is sepa- 
rated in rapidity from the proton, that scatters elastically (see Fig. 1). The t-integrated ep 
cross-section reads: 

da 2na 2 



cm 



[i + (i-y) 2 ]*? i3) ((3,Q 2 ,t) (i) 



d/3dQ 2 d£ /3Q 4 

where the invariants read y = (kq)/(kP), Q 2 = -q 2 , £ = (Q 2 + M X )/{W 2 + Q 2 ) and 
t = (P' — P) 2 . The quantity W 2 = (P + q) 2 is the invariant mass squared in photon-proton 
scattering, and M\ is the invariant mass of the hadronic state X. The reduced-cross-section 
may be expressed in terms of the diffractive structure functions 

o^C&tfU) = F? i3) + J~!^X (3) > (2) 

whereas the structure functions T, L may be, respectively, expressed through transversally 
and longitudinally polarized 7* - proton cross sections = (Q A /4:iT 2 a em P^)da1j'/dM x . 

In the recent analysis jij] the ZEUS diffractive data were fitted within NLO DGLAP 
approximation. A satisfactory description was found only for Q 2 > Q min = 5 GeV 2 . The 
ZEUS fits were performed above <3 mm and then extrapolated to lower photon virtualities. 
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Fig. 2. Left panel - the quark box contribution. Right panel - the qqg contribution. 

The deviations of the fits rapidly grow with decreasing £ and Q 2 reaching 100 percent effect 
at the minimal Q 2 = 2.5 GeV 2 and £ ~ 4 • 10~ 4 . We confirmed this result throug the 
calculation of x 2 /d.o.f. for subsets of ZEUS LRG data with Q 2 > Q 2 nin and (3 > 0.035 
(see Fig. 1, right panel). The cut-off in /3 is imposed to reject part of the data with 
significant contributions from higher Fock states not included in our model. It is clear from 
this discussion that the leading twist DGLAP evolution is unable to describe the DDIS data 
below Q 2 ~ 5 GeV 2 and at the low £. 



III. ESTIMATION OF THE HIGHER TWIST CONTRIBUTIONS 

The large energy limit of the DDIS scattering may be described within the framework 
of the colour dipole model Q, 0]. In this approach the 7*p process is factorized into an 
amplitude of photon fluctuation into the partonic debris and then scattering of these states 
off the proton by the multiple gluon exchange. We take into account the contributions 
from the fluctuation of the photon into a colour singlet quark-antiquark pair qq and into 
gg-gluon triple (see Fig. 2). This gives the ^-integrated 7*p cross section do\^jdM\ = 
daf T jdM 2 x + daf g T /dM 2 x . 

Assuming an exponential t-dependence of diffractive cross-section, one finds for the qq 
component (see Fig.2, left panel) 



da 



qq 

L.T 



fi re>^,;,JQ,z,F)a d (r,t) 



hh,X 



(3) 



where bo is a diffractive slope, zz = z(l — z) and the first sum runs over the three light 
flavours. The second sum of ([3]) means summation over massless (anti)quark helicities (h)h 
in the case of longitudinal photons whereas for transverse photons there is an additional 
average over initial photon polarizations A. The squared photon wave functions can be 
found in literature [f| . 



We use the GBW parametrization |5j for the dipole-proton cross section a^r, £) = 
cr (l — exp(— r 2 /4i?|)) where the saturation radius in DDIS = (£,/x ) x / 2 GeV -1 and 
o"o = 23.03 mb, A = 0.288, xq = 3.04 ■ 10~ 4 . The contribution of the qqg component of 7* 
(see Fig. 2, the right panel) is calculated at (3 = and in the soft gluon approximation 
(the longitudinal momentum carried by the gluon is much lower then carried by the qq 
pair). This approximation is valid in the crucial region of M\ Q 2 or <C 1, where the 
deviations from DGLAP are observed. The correct /9-dependence is then restored using a 
method described by Marquet [7], with kinematically accurate calculations of Wsthoff [8]. 
With these approximations one obtains: 
da™ 1 N c a s a 2 



dMl 167r6n 2tt 2 Ml 
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N 2 qqg (r i) = I d 2 ^-^- (N 02 + N 12 - N 02 N 12 - N ( 



01; 



02' 12 



where Nij = N(fj —fty, roi,ro 2 ,ri 2 = r 02 — r m denote the relative positions of quark and 
antiquark (01), quark and gluon (02) in the transverse plain. The form of N 2 follows from 

n 

the Good- Walker picture of the diffractive dissociation of the photon (9j. The factor \jM\ 
is a remnant of the phase space integration under the soft gluon assumption. The twist 
decomposition of ([3]) is performed through the Taylor expansion in the inverse powers of 
QR whereas that of (jl]) using Mellin transform technic 



IV. DISCUSSION 



In Fig. 3 we compare selected results with data. The saturation model (MSS model) 
results are obtained using the original GBW parameters A and <To, and three massless quark 
flavours. In our approach we modified the GBW parameter xq to £0 = 2xo in order to 
account for the difference between Bjorken x and pomeron £, the variables used in GBW 
dipole cross-section in DIS and DDIS respectively. We chose ot s = 0.4 that provides a 
good description of the data. The conclusion from the analysis and from Fig. 3 is that a 
combination of the DGLAP fit and twist-4 and twist-6 components of the model gives a 
good description of the data at low Q 2 . Inclusion of these higher twist terms improves the 
fit quality in the low Q 2 region (see the dashed curve at Fig. 1 right panel). Indeed, the 
maximal value of x 2 /d.o.f. ~ 1.5 at Q 2 min = 2 GeV 2 is significantly lower then x 2 /d.o.f. ~ 3 
of the DGLAP fit. Nevertheless, it is important to stress that a truncation of the twist series 
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Fig. 3. The LRG ZEUS data for £o> D(3) at low Q 2 compared to a DGLAP fit [2] and the DGLAP 
fit with included twist-4 and twist-4 and 6 corrections from the MSS saturation model. In yellow 
(gray) — the region of (3 where the correction due to qqgg may be neglected. 

(up to twist-6) is required to have a good description of the data. The truncation of this kind, 
however, may be motivated in QCD. Let us recall that in BFKL, at the leading logarithmic 
approximation, only one reggeized gluon may couple to a fundamental colour line. Since 
DGLAP and BFKL approximations have the same double logarithmic (lnxlnQ 2 ) limit, one 
concludes that also in DGLAP couplings of more than two gluons to a colour dipole is much 
weaker than in the eikonal picture. Thus one can couple only two gluons to a colour dipole 
and up to four gluons to qqg component (two colour dipoles in the large N c limit) without 
BFKL constraint. This means that one may expect a suppression beyond twist-8 if only the 
qq and qqg components are included in the calculations. 

In conclusion, the DDIS data at low Q 2 provide the first evidence for higher twists 
effects in DIS in the perturbative domain and opens a possibility for further theoretical and 
experimental investigations. 

Acknowledgments The work is supported by the Polish National Science Centre grant 
no. DEC-2011/01/B/ST2/03643. 

References 



[1] J. Collins, Phys. Rev. D57 (1998) 3051. 



6 



[2] S. Chekanov et al. [ZEUS Collaboration], Nucl. Phys. B831 (2010) 1; Nucl. Phys. B816 (2009) 
1. 

[3] L. Motyka, M. Sadzikowski and W. Slominski, arXiv:hep-ph/1203.5461 (2012). 
[4] N. N. Nikolaev and B. G. Zakharov, Z. Phys. C49 (1991) 607. 

[5] K. J. Golec-Biernat and M. Wusthoff, Phys. Rev. D59 (1998) 014017, Phys. Rev. D60 (1999) 
114023. 

[6] H. Kowalski, L. Motyka and G. Watt, Phys. Rev. D74 (2006) 074016. 

[7] C. Marquet, Phys. Rev. D76 (2007) 094017. 

[8] M. Wusthoff, Phys. Rev. D56 (1997) 4311. 

[9] S. Munier and A. Shoshi, Phys. Rev. D69 (2004) 074022. 



